Majorization theory approach to the Gaussian channel minimum entropy conjecture 
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A longstanding open problem in quantum information theory is to find the classical capacity of 
an optical communication link, modeled as a Gaussian bosonic channel. It has been conjectured 
that this capacity is achieved by a random coding of coherent states using an isotropic Gaussian 
distribution in phase space. We show that proving a Gaussian minimum entropy conjecture for a 
quantum-limited amplifier is actually sufficient to confirm this capacity conjecture, and we provide 
a strong argument towards this proof by exploiting a connection between quantum entanglement 
and majorization theory. 
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During the 1940s, Shannon developed a mathematical 
theory of the ultimate limits on achievable data trans- 
mission rates over a communication channel a work 
that has been central to the advent of our information 
era. Since information is necessarily encoded in a physi- 
cal system and since quantum mechanics is currently our 
best theory of the physical world, it is natural to seek 
the ultimate limits on communication set by quantum 
mechanics. Since the 1970s, scientists started investi- 
gating the improvements that quantum technologies may 
bring to optical communication systems, see e.g. 0t1|- 
Because no proper quantum generalization of Shannon's 
theory existed at that time, the usual approach was to 
compare the performance of different encoding and de- 
coding schemes for a given optical channel. This provides 
lower bounds, but does not give the ultimate capacity nor 
the optimal quantum encoding and decoding techniques. 

In the 1990s, Holevo, Schumacher and Westmoreland 
[H> [H > set the basis for a quantum generalization of Shan- 
non's communication theory. Consider a quantum chan- 
nel M. and a source A = {p a , p a } of independent and 
identically distributed (i.i.d.) symbols. For each use of 
the channel M., Alice sends the quantum state p a with 
probability p ai encoding the letter a. One defines the 
Holevo information 



X (A, M) = S (M(pj) - Y.P« S ( M M) 



(1) 



where p = ^}2 a PaPa and S(p) is the von Neumann en- 
tropy of the quantum state p [7J] . The Holevo information 
X gives the highest achievable communication rate over 
the channel Ai for a fixed source A, which may require 
a collective quantum measurement over multiple uses of 
the channel in order to achieve the optimal decoding op- 
eration. By maximizing Eq. (fTJ) over the ensemble of 
i.i.d. sources A under an energy constraint, we obtain 
the Holevo capacity 



For some highly symmetric channels, such as the qubit 
depolarizing channel, the Holevo capacity actually gives 
the ultimate channel capacity. For a long time, it was 
widely believed that this situation prevails for all chan- 
nels, that is, it was assumed that input entanglement 
could not improve the classical communication rate over 
a quantum channel. However, this was disproved in 
so that the best definition of the classical capacity that 
we currently have requires the regularization 



C(M) = lim -C H (M® n ). 

n— >oo fi 



(3) 



C H (M) = maxx(AM). 



(2) 



where A4® n stands for n uses of the channel. 

An important step towards the elucidation of the clas- 
sical capacity of an optical quantum channel was made 
in [9(, where the authors showed that C(A4) of a pure- 
loss channel — a good (but idealized) approximation of an 
optical fiber — is achieved by a single-use random coding 
of coherent states using an isotropic Gaussian distribu- 
tion. It had long been conjectured that such an encod- 
ing achieves C(M) of the whole class of optical chan- 
nels called phase-insensitive Gaussian bosonic channels 
9], including noisy optical fibers and amplifiers. Actu- 
ally, proving a slightly stronger result known as the min- 
imum output entropy conjecture, namely that coherent 
states minimize the output entropy of phase-insensitive 
channels, would be sufficient to prove this conjecture on 
the capacity of such channels [l(|. Unfortunately, both 
conjectures have escaped a proof for all phase-insensitive 
channels but the pure-loss one. 

In this Letter, we attempt to prove the minimum out- 
put entropy conjecture for a single use of an arbitrary 
phase- insensitive Gaussian bosonic channel Ai, which is 
believed to capture the hard part of the conjecture for 
multiple uses of the channel. We show, using a decom- 
position of any phase-insensitive channel into a pure-loss 
channel and a quantum-limited amplifier, that solving 
the conjecture for a quantum-limited amplifier is suffi- 
cient. This opens a novel way of attacking the conjec- 
ture, using the Stinespring representation of an amplifier 
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FIG. 1: Any phase-insensitive Gaussian bosonic channel M 
is indistinguishable from a composed channel A o £., where 
C is a pure-loss channel and A a quantum-limited amplifier. 
The Stinespring dilation of £ is a beam-splitter of transmis- 
sivity T, while the amplifier A of gain G becomes a two-mode 
squeezer of parameter r (G — cosh 2 r) in which the input 
mode A interacts with a vacuum environmental mode E. 



channel as a two-mode squeezer, and exploiting the con- 
nection between entanglement and majorization theory. 

Quantum model of optical channels.- A quantum opti- 
cal channel can be modeled as a Gaussian bosonic chan- 
nel. It is a trace-preserving completely positive map fully 
characterized by the action on the Weyl operators of two 
2x2 matrices, K and N [ll| - [l3j . An intuitive under- 
standing of K and N is given by the action of the channel 
on the mean x and second moments 7 of the input state, 



x — > Kx, 7 — > KjK 



N, 



(4) 



For the map to be completely positive, K and N must 
satisfy [lj] 



N>0, detiV > (detK - l) 2 



(5) 



resulting from spontaneous emission during the amplifi- 
cation process; (iii) the additive classical noise channel, 
corresponding to r = TG = 1 and added thermal noise 
n = 2G - 1. 

Reduction of the minimum entropy conjecture.- As 
stated earlier, our ultimate goal is to address the fol- 
lowing conjecture: 

CI. Coherent input states minimize the out- 
put entropy of any phase- insensitive Gaussian 
bosonic channel M.. 

Three simplifications can be made at this point. First, 
due to the concavity of the von Neumann entropy, the 
minimization can be reduced to the set of pure input 
states. Secondly, applying a displacement D(a) at the 
input of the channel has the same effect as applying 
D(- s /ra) at the output, i.e., MoD(a) = D(y/ra)oM. So, 
because the von Neumann entropy is invariant under uni- 
tary evolution, we can restrict our search to zero-mean 
input states, that is, states \<p) satisfying (<£>|a|<^) = 
where a is the modal annihilation operator. Finally, ex- 
ploiting the decomposition M — A o C it is easy to see, 
using the concavity of the von Neumann entropy, that the 
minimum output entropy of channel M is lower-bounded 
by that of channel .A, i.e., min^S(M((j))) > min^S(A(ip)) 
[la ]. Since the vacuum state is invariant under £, we con- 
clude that proving that vacuum minimizes the output en- 
tropy of channel A implies that vacuum also minimizes 
the output entropy of channel M. 

The previous straightforward derivation shows that the 
conjecture CI is strictly equivalent to the following one: 

C2. Among all zero-mean pure input states, 
the vacuum state minimizes the output en- 
tropy of the quantum-limited amplifier A. 

Entanglement and majorization theory.- The Stine- 
spring dilation of a quantum-limited amplifier of gain G 
is a two-mode squeezer of parameter r, with G — cosh 2 r, 
which effects the unitary transformation (see Fig. [TJ 



Most naturally occurring optical channels, such as op- 
tical fibers or amplifiers, are phase insensitive. They 
correspond to K = diag(- v /r, y/r) and N = diag(n, n), 
with t being either the attenuation < r < 1 or the 
amplification 1 < t of the channel, and n being the 
added noise variance. Using the composition rule of 
Gaussian bosonic channels [lj], it is easy to show that 
every phase-insensitive channel M. is indistinguishable 
from the concatenation of a pure-loss channel C of trans- 
missivity T with a quantum-limited amplifier A of gain 
G, see Fig. Q] The parameters T and G must satisfy the 
relations r = TG and n = G(l - T) + (G - 1) in order 
to guarantee M. = A o C . Three limiting cases are of 
particular interest: (i) the pure-loss channel, correspond- 
ing to G = 1 and < T < 1, having a quantum- limited 
noise of n = 1 — T; (ii) the quantum-limited amplifier 
corresponding to T = 1 and G > 1 , with noise n = G — 1 



U{r) ~ exp r(aAdE 



a\a E 



L)/2 



(6) 



between the input mode A and an environmental mode 
E, where a) x and ax are the creation and annihilation op- 
erators of mode X. Because the entanglement E [IV^ab] 
of a pure bipartite state \ip) ae is uniquely quantified by 
the von Neumann entropy of its reduced density opera- 
tor p A = Tr B [|^) AB (0|], i.e., E [\4>)ae] = S (pa), we can 
equivalently rephrase conjecture C2 as 

C3. Among all input states \4>)ae = \<f) <8> |0) 
of a two-mode squeezer with \ip) having a zero 
mean, the vacuum state \0}ae = |0) <g> |0) 
minimizes the output entanglement. 

In the remainder of this Letter, we exploit the con- 
nection between entanglement and majorization theory 
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to attack the proof of C3. Majorization theory pro- 
vides a partial order relation between probability dis- 
tributions [l5|, HB). One says that a probability distri- 
bution p = (pi, ...,pd) T (d might be infinite) majorizes 
another one q (denoted p >- q) if and only if there exists 
a column-stochastic matrix D (a square matrix whose 
columns sum to one) such that q = Dp, showing that 
q is more disordered than p. It implies that all concave 
functions of a distribution, most notably the entropy, can 
only increase along such a "disorder enhancing" transfor- 
mation. From an operational point of view, an interesting 
way of proving majorization is by checking the relations 

m m 

XX>XX V ™< d > (7) 
n—1 n—1 

where and q^ are the original vectors with their com- 
ponents rearranged in decreasing order. The notion of 
majorization can be extended to entangled states [17| : a 
bipartite pure state \<j>) majorizes another one (noted 
\(j)) >~ \ifj)) if and only if the Schmidt coefficients of \<j>) 
majorize those of This guarantees the existence 

of a deterministic protocol involving only "local opera- 
tions and classical communication" (LOCC) that maps 
\tp) into ensuring the relation S[|^/))] > 25[|0)]. We 
are now ready to introduce the following stronger conjec- 
ture (it implies C3): 

C4. For any zero-mean state \ip) ^ |0), the 
state U(r)(\(p) <g> |0)) is majorized by the two- 
mode squeezed vacuum state J7(r)(|0) €5 10)). 

Infinitesimal two-mode squeezer.- Before addressing 
the general case, let us prove C4 for an infinitesimal two- 
mode squeezer by expanding the unitary transformation 
(|6|) to the first order in the squeezing parameter r, 

U{r)=I+ r -(a A a E -a\a^+0{r 2 ), (8) 

where I is the identity operator. Defining the state 
|</j_l_) = —a' A \<p)/(l + n^) 1 / 2 , where n v = ((p\a A a,A\(p) 
is the mean photon number of the input state \ip), the 
output state becomes 

\4>out) AE ~ 

y/\>W) ® 1°) + v 71 - ® |i), (9) 

with X v = 1/[1 + r 2 {n v + l)/4]. For any physical state 
\ip) with finite energy n v , one can choose r small enough 

_ 1/2 

so that the condition r <^C 1 is satisfied and the ap- 
proximation @ holds. The key point is to realize that 
since the input state \(p) has a zero mean, the states |<p_i_) 
and | ip) are orthogonal, so that the state ((9|) is already 
in Schmidt form. Therefore, if \ip) and \tt) are two input 
states such that n v < fi n , then X v > A w , implying that 
U(r)(\<p)®\0)) y U(r)(\n)®\0)) as a result of Eq. . In 
other words, any output state is majorized by the states 
having a lower mean input photon number. Finally, since 
the vacuum state has the minimum mean photon number 



(fi v = 0), this majorization relation proves conjecture C4 
for infinitesimal two-mode squeezers. 

Majorization relations in a two-mode squeezer.- In or- 
der to address the conjecture C4 for any r, let us consider 
the number-state expansion of an arbitrary input state 
\ip) = YlkLo c k\k), which leads to the output state 

oo 

CT(r)(b)®|0))=£c fc |*W), (10) 
k=a 

where A = tanhr and |v?^ ) stands for the output state 
corresponding to an input Fock state | ip) =\k). As shown 
in [l5j], we have 



^ I 

K fe) > = E ^Jp [ n\X)\n + k)®\n), (11) 

n=0 

with Schmidt coefficients 

p W(A) = (l-A 2 ) fe+1 A 2 «^ + fc ). (12) 

We have been able to prove two chains of majorization 
relations by considering either different Fock states \k) at 
the input (for a fixed squeezing parameter r) or different 
values of r (for a fixed input Fock state First, when 
restricting to Fock states \k), we can prove that 

y (13) 

since there exists a column-stochastic matrix 

D nm = (1 - A 2 )A 2 <"-" l )ff(n - m), (14) 

such that p( fc+1 ) (A) = DpW (A), where H(x) is the Heav- 
iside step function defined as H(x) = for x < and 
H(x) = 1 for x > 0. The details of the proof are pro- 
vided in (l5j . where we also give the explicit form of an 
LOCC protocol that deterministically maps into 
1*-^ )• Iterating this procedure, we can easily prove that 

I^A^) ^~ )■> > k, for which we also give the cor- 
responding column-stochastic matrix and deterministic 
LOCC protocol. 

For our matters here, the central consequence is that 

l*i 0) ) >~ I*!?): ^ k > °) that is > we nave proved con- 
jecture C4 for the restricted, but complete, set of input 
Fock states. Remarkably, this would be sufficient to prove 
the single-use minimum entropy conjecture if it could be 
shown that the output-entropy minimizing input state 
is isotropic, i.e., its Wigner distribution is rotationally 
invariant. This is because the Fock states are the only 
isotropic, zero-mean pure states. 

Second, given an input Fock state one can show 
that there exists a majorization relation in the direction 
of decreasing squeezing parameter, that is, 

|*W)H*i fe) } VA'<A, (15) 
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FIG. 2: Entanglement of the output state as a function 

of the squeezing parameter r. As explained in the text the 
entanglement is monotonically increasing with r for all Fock 
input states, while, for a fixed r, it monotonically increases 
with k. This behavior is in full agreement with the majoriza- 
tion relations (|13|l and (|15p proved in the text. The arrows in 
the figure indicate the majorization order. 



since one can build 1 1511 a column-stochastic matrix 



R (k) = = 



m + k 
m 



T (k,m) ,2 



1-A 2 
1 - A' 2 



H(n - m) (16) 



- (fe,m+l) w2 
J n-m-l A 



A 2(r 



with 



of the first 21 Fock states, which strongly suggests that 
conjecture C4 holds. 

Conclusion.- Using the decomposition of phase- 
insensitive Gaussian bosonic channels into a pure-loss 
channel and a quantum-limited amplifier, we have shown 
that proving a reduced conjecture for the quantum- 
limited amplifier is sufficient to prove the single-use min- 
imum entropy conjecture. Using Stinespring's theorem, 
this boils down to proving that the vacuum minimizes the 
output entanglement of a two-mode squeezer. Then, us- 
ing the connection between entanglement and majoriza- 
tion theory, we have provided a partial proof of this con- 
jecture for a special class of input states, namely photon 
number states, as well as a full solution for the infinites- 
imal channel. To prove the conjecture in general, we 
are left with the (possibly simpler) task of showing that 
the output-entropy minimizing input state is isotropic in 
phase space, that is, no symmetry breaking occurs. Thus, 
apart from reinforcing the conjecture even further, we be- 
lieve that our analysis offers a new possible approach to 
its proof. 
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X'-^B(X u ;n,l + k), (17) 
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and B(z;a,b) — J* dx x a 1 (1 — x) b 1 being 
the incomplete beta function, such that pw(A) = 
R {k) (\, A')p (fe) (Y)- In [H| we give a deterministic LOCC 



(*0 



I* 



0)\ 

A' 



protocol performing the transformation |W 

In Figl21 we summarize the two chains of majorization 
relations and their implications on the output entangle- 
ment. From this, as well as the case of the infinitesi- 
mal two-mode squeezer, it is tempting to conclude that 
n v <n- K always implies U(r)(\(p)<g)\0)) y U(r)(\n) ® |0)). 
However, we have numerically observed that this does 
not hold in general, which probably reflects the diffi- 
culty of proving the conjecture. As a concrete example, 
we note that the state U (r)[(y/0A\l) + VoH\2)) ® |0)] 
has n = 1.6 mean input photons but is less entangled 



for r > 0.75 than I* 



Nevertheless, our numerical 



investigations have shown that for an arbitrary input 
state | ip), the output states corresponding to different 
squeezing parameters satisfy the majorization relation 
U(r'){\<p) <g) |0)) y U(r)(\tp) <g) |0)) for r' < r. Further- 
more, we have numerically checked that for a fixed r, the 
majorization relation U(r)(\0) <g) |0)) y U(r)(\<p) <g> |0)) is 
satisfied by tens of thousands of random superpositions 



In what follows, we give a more complete overview of 
the calculations leading to the main results of this Let- 
ter. First, we derive the lower bound used to reduce 
conjecture CI to C2. Second, we review the concept of 
majorization in probability theory, and describe its use 
in the context of quantum entanglement. Then, we de- 
tail the calculation of the output state of a two-mode 
squeezer for an arbitrary input state expressed as a su- 
perposition of Fock states. Finally, we provide a detailed 
derivation of the chain of majorization relations that are 
obeyed by a two-mode squeezer with number-state inputs 
in one port, and present their associated local operation 
and classical communication (LOCC) protocols. 



1. Reduction of the minimum entropy conjecture 

In what follows we exploit the decomposition A4 = 
A a C and the concavity of the von Neumann entropy 
to prove that the minimum output entropy of chan- 
nel A4 is lower-bounded by that of channel A, i.e., 
min S'(.M (</>)) > rwD4,S(A(ip)). 
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Let \<p) be an input pure state of channel A4. After 
passage through the pure-loss channel £, the inter midi- 
ate state (between C and A) is a = C(\4>) {<p\) ■ For any 
decomposition {pi,^} of a satisfying a = j^Pil^iX^il) 
we have the following chain of inequalities 

S(M(\<f>)(4,\)) ( = } S(A(a))®sfe Pi Amm) 

(3) 

i 

(4) 

> mm.^S(A(^)), (A.l) 

where we have used: the channel decomposition M = 
A o L in (1); the linearity of quantum operations in (2), 
the sub-additivity of von Neumann entropy in (3); and, 
finally, the definition of the minimum output entropy of 
channel A in (4). The proof concludes by noticing that 
Eq. (|A.1[) holds for every input state of channel A4, in- 
cluding the one minimizing the output entropy of M. 

2. Majorization and Entanglement 

Majorization appeared as a way to order probability 
distributions in terms of their disorder, in an effort to 
understand when one distribution can be built from an- 
other by randomizing the later [la ] . Take two probability 
vectors p = (pi,p2, —,Pd) T and q = (qi,q 2 , -,qd) T of di- 
mension d (which can be infinite as in our case), properly 

normalized, that is, Yln=iPn = Er>=i In = L We say 
that p majorizes q, and denote it by p >- q, if and only 
if 

m m 

$>!l>£«£ Vm<d, (A.2) 

n—l n—1 

where p^ and q^ are the original vectors with their com- 
ponents rearranged in decreasing order. This definition 
is useful from a practical point of view, since it is easy to 
check numerically if two vectors satisfy these relations. 
Nevertheless, it can be proven that p >- q is strictly 
equivalent to two other operational relations: 

Ml. For every concave function h{x), we 

have En=l h (Pn) < En=l Kin)- 

M2. q can be obtained from p via q = Dp, 
where D is a column-stochastic matrix. 

A square matrix D is column-stochastic if its elements 
are real and positive, its columns sum to one, and its 
rows sum to less than one. Most of the literature on the 
connection between majorization and quantum informa- 
tion studies finite-dimensional systems, in which case it 
can be shown that column-stochastic matrices are also 
doubly-stochastic (columns and rows both sum to one). 



In this work we need the slightly more general definition 
of column-stochastic to cope with infinite dimensional 
spaces (l8| . Physically, stochastic matrices are equiva- 
lent to convex mixtures of permutations of the vector 
components, and hence, property M2 shows that q is 
more disordered than p. 

Interestingly, majorization theory can also be used to 
answer the question of whether Alice an Bob can trans- 
form a shared bipartite pure state \iP)ab into \tp)AB by 
using a deterministic protocol involving only local op- 
erations and classical communication (LOCC) [l7], Il9j - 
Given the probability vectors p^, and p v generated with 
the Schmidt coefficients of these states (the eigenvalues 
of the reduced density operators) , it is possible to prove 
that the transformation \4>}ab — ► \*-p) ab is possible if and 
only if p v >- p^, that is, if the Schmidt coefficients of 
\tp) ab majorize those of \iP)ab, in which case we use the 
symbolic notation \(p)ab >~ IV") ab- The entanglement 
of a pure bipartite state \iP)ab being measured by the 
von Neumann entropy of the reduced density operator 
Pa = Tr^l^As], and the von Neumann entropy being a 
concave function, one gets as an intuitive corollary that 
\iP)ab can only be transformed deterministically by an 
LOCC protocol into states of lower entanglement, i.e., 

EU)ab] > E[\<p) AB ], (A.3) 

as follows from property Ml. 

Note that while \<p)ab >~ \iP)ab implies that p v can 
be transformed into p^, by application of a column- 
stochastic matrix, the transformation goes in the oppo- 
site direction for the corresponding states, that is, it is 
\ip) A B the state which can be transformed into \^p)ab by 
a deterministic LOCC protocol. In other words, at the 
level of probability distributions the transformation in- 
duces disorder (increases the entropy), while at the level 
of states the transformation decreases the entanglement, 
as corresponds to physical deterministic LOCC protocols. 

3. Output States of a Two-Mode Squeezer 

If we inject the vacuum state at the input of a two- 
mode squeezer U(r), we obtain the two- mode squeezed 
vacuum state 

|*(°)) = U(r)\0,Q) = — — VtanW \n,n), (A.4) 
coshr ^-^ 

where \n) is a number state, and we use the compact 
notation \rn)A ® \n) B = \m, n). 

Consider now the more general input state 

oo 

\4>) = \ v )®\Q) = Y J c n \n,Q), (A.5) 

n=0 

written in the number state basis, which becomes the 
state 

oo 

\4>out)=U{rm = Y^c n \^ {n) ), (A.6) 

n=0 



6 



with 



with 



|*W) = tf(r)|fc,0>, 



(A.7) 



after passing through the two-mode squeezer. 

In the reminder of this section, we focus on finding 
a manageable expression for the states I^W), that is, 
for the output state of the two-mode squeezer when a 
number state \k) is fed through one of its input ports. 
We start by noting that |\I>( fc )) can be written in terms 
of the two- mode squeezed vacuum state l^ )) as follows 

|*<*>) = J=tf(r)a>,0) = -j=lU(r)a\u(ry} k \¥V), 

(A.8) 

which, using the relation 

U{r)a\U{r)^ = coshr a) A — sinhr as, (A. 9) 

can be rewritten as 



pW(A) = (l-A 2 ) fc + 1 A 2 " /5/ W '' 



(A.17) 



to stress the dependence of the state on the squeezing pa- 
rameter. Note that the states (|A.16|) are already written 
in Schmidt form, and in the following we will use 



p<*> = (p?U fc) ,...) T , 

to denote the corresponding probability vectors. 



(A.18) 



4. Proof of the Majorization Relations for Fock 
State Inputs 

In this section we will explain how to derive the 
column-stochastic matrices needed to prove the majoriza- 
tion relations employed in the Letter. 



|tf<*)) = J2 ( ^ ' ( fc ) coshV sinner a^a k B ] \^). 
j=o \JJ 

(A.10) 

Now, an easy calculation shows that 

1 oo 

,|*(°)) = — — Vv^itanhV \n,n- 1) (A.ll) 

pn«n r — ' 



as I 



= — y/m~+T tanh m+1 r |m + 1, m), 

n->m+i coshr ^ — ' 



leading to the following identity 

a s |* (0) ) = tanhr a\\¥ 0) ), 
which allows us to rewrite (|A.10[) as 



(A.12) 



^( fc)) = ^Lj2(-l) k - j ( k ) tanh 2 ^Va^'|*(°)). 

(A.13) 

Finally, using the relations 

£(-l) fc -^V- j - (1-*)*, (A.14a) 

3=0 

l-tanh 2 r = cosh -2 r, (A. 14b) 
we can write the previous expression as 
1 



/fc!cosh fc r A 



cosh fe+1 r5; 



tfc|^(0)v 



n + k 



(A.15) 



tanh"r |n + k,n). 



Let us define A = tanhr; from now on we will use the 
notation 



*! fc) > = Evp« ) ( A )i™+^>' 



(A.16) 



Proof of 



Because the states ) are already in Schmidt form 
as commented previously, we need to prove that there 
exists a column-stochastic matrix D such that 



(k) 



(A.19) 



This is actually quite simple if one notices that the Pascal 
identity 



n + k + 1 
fc + 1 



n + k 
k 



n + k 
k + l 



(A.20) 



implies the following relation (with the convention pffi 
for n < 0): 



(i-A 2 )p«+A 2 P ri 



2 



(A.21) 



This recurrence allows us to connect p( fc+1 ) with p( fe ) by 
means of a lower-triangular matrix 



p[ k+1) 



(1-A 2 ) 



/ 1 

A 2 1 

A 4 A 2 1 

V ; 



P [ k) 
pf ] 



I V 



/ 

(A.22) 



(A.23) 



n=0 



or in a more compact notation 

n 

P r = E(i-A 2 )A 2 >ii. 

m— 

It is fairly easy to show that the triangular matrix shown 
above, whose elements are explicitly given by 

D nm = (1 - A 2 )A 2 ("- m ^( ? i - m), (A.24) 

with H(x) being the Heaviside step function defined as 
H(x) — 1 for x > and H(x) = for x < 0, is column- 
stochastic. Hence we conclude that >~ |^ fc+1 ^) as 
commented in the Letter. 
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b. Proof of y |*( fc+Afc) ) for Ak > 



I* 



(fe) 



(fe+1) 



implies | ^ A 



I* 



(fc)\ 



>- 



l^^-rii/v;^ a ^ ^ o (note that Afc is a positive inte- 



It is clear that 

,0+Afc) 

ger by definition), as majorization is clearly a transitive 
relation. This shows that when restricted to Fock-state 
inputs, the output entanglement of a two-mode squeezer 
increases monotonically with the number of input pho- 
tons. 

In order to find the explicit column-stochastic matrix 
jj(Ak) satisfying p( fc + Afc ) — £)( A ' £ )p( fe ) i we use the inde- 
pendence on fc of the matrix D which allows us write 



D 



(Afc) _ 



D x D x ... x D. 



(A.25) 



Afc times 



An explicit form of the elements of this matrix can be 
inferred for any Afc by evaluating the first matrices: 



D 



(2) _ 



2 x 2 



(1 - A ) 



D 



(3) _ 



2\3 



= (1-A 2 ) 4 



/ 1 0. 

2A 2 1 0. 
3A 4 2A 2 1 . 
4A 6 3A 2 2A 2 1 . 

V ; ; i ; • 

( 1 

3A 2 1 

6A 4 3A 2 1 

10A 6 6A 2 3A 2 1 

V ; i ; 

( 1 

4A 2 1 

10A 4 4A 2 1 

20A 6 10A 2 4A 2 1 

V : ! : \ 



squeezing parameter A. In what follows we prove a 
stronger result, that a given two- mode squeezed vacuum 
state majorizes all the two-mode squeezed vacuum states 
with stronger squeezing. 

We seek for a column-stochastic matrix R(X, A') satis- 
fying 



p(°)(A)= J R(A,AV° , (A') 



(0)/ 



(A.28) 



Based on the matrices of the previous sections, we make 
an ansatz in which R is a lower-triangular matrix whose 
columns are all built from a vector r(A, A'), that is, 



R = 



{ r Q 

ri ro 

r 2 n r 

?"3 r 2 n r 

V • ; ; i 



\ 



(A.29) 



Introducing this ansatz into equation (IA.28|1 . and recall- 



ing that Pn \x) = (1 — x 2 )x 2n , we get the following set 
of linear algebraic equations 

(1-A 2 ) = (1-A' 2 )r , (A.30) 
(1-A 2 )A 2 = (l-A' 2 )(A' 2 r +r 1 ), 
(1-A 2 )A 4 = (l-A' 2 )(A' 4 r + A' 2 r 1+ r 2 ), 

which can be solved by recursion leading to the solution 

r„ = [ A " - H ( n !) A ' 2 ] A2( " _1) ' ( A - 31 ) 

which can checked, by induction, to be the solution for a 
general n. Note that Y^=o r « = 1 as expected. 



Hence, all Z?( Afe ) matrices have a similar structure, except 
for the (1 — A 2 ) Afc prefactor, and the numbers accompa- 
nying the powers of A 2 in the columns, which are given 
by the Afcth diagonal of the Pascal triangle. It is then 

fairly simple to prove by induction that the elements of 
£)(Afc) are gj ven by 



= (1 - A 2 )** ( m + fc Afc ; X ) A 2 (— )H(n - m). 



'm + Afc - 1 N 
Afc - 1 

(A.27) 

Note that this general majorization relation implies 
in particular that | V I / A ' ) ) >- ) Vfc, and therefore, 

among all Fock state inputs, the vacuum state is the one 
which minimizes the output entanglement of a two-mode 
squeezer. 



c. Proof of |*™) >- j*™) for A' < A 

It is well known that the entanglement of the two-mode 
squeezed vacuum state monotonically increases with the 



d. Proof of 1*^) >- for A' < A 

The same kind of majorization relation can be proved 
for any |\1/ A ' C ^ ' ) ) state, although the proof is now a little 
more involved, as we need to find a matrix R^(\, A') 
satisfying 



(*0\ 



p (fc) (A) =i? (fc) (A,A')p (fe) (A'), 



(A.32) 



which now depends on the value of fc. As we now prove, 
the matrix RW(\, A') can still be chosen to be lower- 
triangular, but now every column is defined by its own 
vector r( fe 'ft , that is 



' 'o 







Jfe.o) (fc,l) 

'i 'o 

(fc.O) (fc,l) (fc.2 

r, r) r, ' 



r^ k ' ^ r. 



(fc,i) 







„(M) Jk,3) 



'o 



(A.33) 
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Because we have to recover the case k = (|A.31I) . we 
make the following ansatz 



r (k,n) _ A 2(m-1) 



1-A 2 



with B 



(0. 



B 



(k^0,n) 



a 



1-A' 2 

B^X 2 - C^ n) H{m - 1)A' 2 

(o,«; 



(A.34) 



1, and where the coefficients 



and Cm ' may depend on A and A'. 



Similarly to the previous section, we can find the co- 
efficients Bm"' and Cm by introducing this ansatz 
in (| A.32I) . and using the explicit form of the probability 

vectors p ( n\x) = (1 - x 2 ) k+1 x 2n ( n + fe ) . Let us show this 
process step by step. 

The system (IA.32[) can be rewritten in a compact form 

as 



P { n ] W = E rL fe '"- m) (A,A>W ro (A'). (A.35) 



m=0 

For n = 0, this sets 
while for n = lwe get 



R (M) _ -I 

n o — ^ 



(A.36) 



A^j 1 ) - u a- 



2 / " - 1 \ . D (fc,i) */a f fc + ^ _l B^ fi) \ 2 - c[ kfi) X' 2 



1 



(A.37) 

of which = and cf< 0) = B { Q kA) ( k + 1 ) are 

valid solutions. Similarly, for n — 2 (IA.35[) yields 



A 4 f ' + 2 ) = < 2) A' 4 (* + 2 ) + B™ A 2 A' 2 (* + 1 



-C[ k ' l) X' 4 f k | X ] + sf' 0) A 4 - Cf < 0) A 2 A' 2 , (A.38) 



of which B 



(fc,2) 



fc+2\ ^(fe.O) 



( fc + 2 ), a 



c (k,i) = B^ k ' 2) { k + 2 )/( k + 1 ) are now valid solutions. 



B^^f), and 



We observe the pattern of solutions 

B (k,0) = 



^(k.n) _ r>(k,n- 
°m — D m-1 



m + k 
k 

(n+k+l\ 
(fe,n+l) fc ) 



(A.39a) 
(A.39b) 



so that the components of the vectors r^ 8 '™) can be rewrit- 
ten as 



„(fe,n) _ 



fc\ 1 / 1 - A 



1 - A' 2 



(A.40) 



where we have defined the new parameters 



L (k,n) _ (n + k\ {k . 



(A.41) 



which satisfy I/m'°^ = { m ^ k ) except for m < 0, in which 



case L 



(k,- 



0. 



In order to find the coefficients im'"' we use a further 
condition: as R( k \\, A') must be column-stochastic, the 
vectors r' fc '™) must be normalized. Let us then define the 



series 



S (k,n) = L^- n )A' 2m , (A.42) 
m=0 

in terms of which the normalization condition 

m=0 



(fc,n) 



1 can be rewritten as 



y2 g(k,n+l) _ g(k,n) 

Starting from 



n + k 
k 



1-A'' 
1-A 2 



(A.43) 



S (k,0) = g A" + fc A y 2m = (1 _ ^-(fc+l) ^ (A 44) 
m=0 A / 

these relations allow us to find the rest of S^ k,n ' recur- 
sively, obtaining 



S (k,i) = y-2 (1 _ A 2)-(fe+i) [1 _ (1 _ A' 2 ) fc+1 ],(A.45a) 

(A.45b) 



S (k,2) =y -2 (1 _ A 2 ) -(fc+l) | A ,-2 

(1 - A' 2 ) fe+1 



w-2 ik+1 
A + 



S (*,3) = y-2 (1 _ A 2j-(*+i) {y-4 
'fc + 1 \ (k + 2 

* + u 



A'" 4 + A'" 2 



(A.45c) 
(1 - A' 2 ) fc+1 



from which one sees the general pattern 

g(k,n) _ y-2n^ _ A 2-j-(fe+l) 

ri-l 

1 - (1 - A' 2 ) fc+1 E A' 2 ' 

1=0 



(A.46) 



l + k 



The sum on the right-hand side term can be written in 
terms of the incomplete beta function 



B(z;a,b)= dxx a - 1 (l-x) b - 1 , 
Jo 



n-l 

E 

1=0 



l + k 



A' 2i - (1 - A' 2 )-( fc+1 
n + k 



(A.47) 



(A.48) 



1 /•■( " " " )B(X' 2 ;n,k + 1) 



We can therefore rewrite the condition (|A.46|) as 

oo 

E i^' n+1) A 2m (A.49) 

'n + £A 



= A'- 2 "(l-A 2 )- (fc+1) n^ 



£(A' 2 ;n,/c + l), 
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which, given the result (|A.44[) . can be satisfied by choos- 



ing 



1) 



(A.50) 

Note that this expression is valid even for n = 0, as 



lim aB(x] a, b) = 1, 

a— >0 



(A.51) 



when 6 is a positive integer. Introducing this expres- 
sion for the Lm' n ^ coefficients in r^ k ' n ^ (|A.40[) . and this 
into (|A.33[) . we get the column-stochastic matrix R(\, A') 
given in the Letter. Hence, we have been able to find 
a stochastic map connecting pW(A') to p( fc )(A), which 
proves the majorization relation >~ \^\) if A' < A. 



5. LOCC protocols 

For completeness, we now give the LOCC protocols 
corresponding to the previous majorization relations. We 
believe that these could offer an alternative (more phys- 
ical) way of attacking the proof of the conjecture for a 
general input state like (|A.5[) . and hence find it appro- 
priate to explain how to build such protocols. 



Transformation -> |*^ ; ) 



which then yields the desired state regardless of m, 

that is, deterministically. Remark that the shift operator 
is trace preserving in the subspace spanned by {\j + m + 
l)}j=o,i,..i which is the support of (Ia ® B m )\^>^ +1 ^) on 
Alice's side. Notice that one can easily build a shift oper- 
ation that acts on Alice's full Hilbert space by appending 
ancillary qubits. 



b. Transformation j^ fe+Afc) ) -» j*^ fe) ) for Ak > 



Similarly as before but exploiting now (|A.27|i . we en- 
gineer the following POVM on Bob's side 



B m — X/ \ 



Pi 



(fe+Afc) 



\l-m)(l\, 



(A.55) 

which, combined with the conditional shift in Alice's side 



A m =J2\l){l + m + Ak\, 



(A.56) 



1=0 



deterministically transforms the state |^ 

,(k) 



(k+Afc)\ 
A / 



into 



\V\ '). Wheneverfc = 0, we obtain the two-mode vacuum 

squeezed state | V E'^ ' 1 ), which is thus at the end of the ma- 
jorization chain, and its entanglement is minimum when 



compared to all other states [$! 



(fc)\ 



Let us assume that Alice and Bob share the bipartite 
state |^ r A fc+1 ), and want to convert it into }• In- 
spired by the recurrence relation (|A.23|i . we propose the 
following LOCC protocol. Bob starts by performing a 
POVM measurement [2(| described by the measurement 
operators 



B„ 



\ 



(1 - A 2 )A 2m pJ 



(k) 



Pi 



(fe+i) 



-\l-m)(l\. (A.52) 



Using Eq. (|A.23|) . it is easy to verify the condition 
J2m=o BmB m = I. After Bob has completed his local 
measurement, depending on the outcome m of the mea- 
surement, the joint state "collapses" to 



c. Transformation |*^ 0) ) -> |*™} for A' < A 

Constructing an LOCC protocol from the stochas- 
tic matrix R(X, A') (|A.29I) which connects p (0) (A') with 
p(°)(A) is not an easy task. Interestingly, we found a 
simpler deterministic protocol achieving the same result. 
Let us first give a probabilistic scheme performing the 
transformation, which we later make deterministic. 

As shown in Figure Bob mixes his mode B with an 
ancillary mode C on a beam-splitter of transmissivity T. 
The initial state is 

oo 

Wabc = l*i 0) ) ® |0> = ^(A)^A>,n,0), (A.57) 



(I A <X> B m ) |tf 



a X \JPn-m\ n + k + 1 > n - m ) 
n—m 
oo , 

X yPn ] \n + k + m+l,nXA.53) 



Then, after Bob has communicated the outcome m of 
his measurement to Alice, she performs the local shift 
operation 



A ro =5^|I){Z + m+l| > 



(A.54) 



where A/"(A) = (1 — A 2 ) 1 / 2 a normalization factor. After 
passage through the beam-splitter, the joint state be- 
comes 



W)abc = Af{\) £ (TA 2 ) n 



/2 



n 7n=0 



1-T 
T 



m/2 



1/2 



' \ I In, n — m, m). (A.. 

i m ; 



Finally, Bob measures the number of photons reflected 
by the beam-splitter. The outcome of the measurement 
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Mode .4 in 



I* 



(0) 
A 



Mode C 



10) 



ModeB hi 




Photon Counter 



FIG. 3: Probabilistic LOCC protocol achieving the transfor- 
mation |4^ 0) ) -> |*™) for A' < A. Initially, Alice and Bob 

share the entangled state \^ ° )ab. The first step of the pro- 
tocol consists in Bob mixing his mode B with a vacuum an- 
cillary mode C into a beam-splitter of transmissivity T, and 
measuring the number of photons at the output of mode C 
with a photon counter. Conditioned to the measurement of 
zero reflected photons, the desired transformation is achieved 
with A' = VTX. 



will be zero with probability V = N 2 (VT\)/Af 2 (\), after 
which the state will collapse according to 



T>W) 



AB 



ABC 



N{\)Y T"/ 2 A"|n,n) 



n=0 



P\9 



(0) 



■^ab. (A.59) 

Then, by choosing the transmissivity of the beam-splitter 
to satisfy A' = VT\ we obtain the target state. Note that 
there always exists a valid transmissivity T, as A' < A. 
The input state \^^}ab ® |0)c being a Gaussian state 
and the projection into vacuum being a Gaussian op- 
eration, there must exist a deterministic LOCC proto- 
col generating the same outcome [ll[. Such a protocol 
consists of replacing Bob's projection onto vacuum by 
heterodyne detection followed by local displacements on 
Alice and Bob sides that are proportional to the outcome 
of Bob's heterodyne measurement. 



Transformation l^?') 



for A' < A 



Similarly to the case k = 0, constructing an LOCC pro- 
tocol from the stochastic matrix i?( fe )(A, A') ([A. 331) which 
connects p( fe )(A') with p( fc )(A) is not an easy task. In- 
stead, we give a simpler deterministic protocol achieving 
the same result. 

Just as in the previous protocol, Bob starts by mixing 
mode B with an ancillary mode C on a beam-splitter of 



transmissivity T. The joint initial state is 



I ABC 



= K fc) >®|0) 

oc 

= AA(M)E AT ' 



(A.60) 



n=0 



1/2 



1 7i + k, n, 0), 



with Af(k, A) = (1 - A 2 )( fe+1 )/ 2 , which becomes 



W)abc=M(k\) F#) n 



/2 



n.m— 



l-T 
T 



i/2 



I ! I J |n + k, n — to, to), (A. 61) 



A- 



after passing through the beam-splitter. 

Second, Bob measures the number of photons reflected 
by the beam-splitter. With probability 



m = {1 _ T y x xf k + l \ -^ 2 (M) 



i J N 2 (k + i,VTxy 



(A.62) 



the outcome of the measurement will be I photons, and 
the state of modes A and B will collapse in that case to 



VW) W')ab = c{iW) 



■■ Af(k, A) 



1 - T 
T 



ABC 
1/2 



Y,(T\ 2 ) n/2 



n—l 



. , \ 1/2 / x 1/2 



(A.63) 



1 7i + k, n — I). 



Now, making the variable change n — I — > n in the sum, 
and using the relation 



7i + I + k\ (n + l\ fn + k + f\ fk + l 



k J \ I J \ a 
this state can be rewritten as 

VWW)ab = M(k,X)(l-T) l ^X l 



(A.64) 



1/2 



^ (TA2yi/2( 'n-^/y /2 |n Lli) 



(A.65) 



n=0 



Notice that by properly choosing the transmissivity of 
the beam-splitter so that A' = v^A, the final state is 
Therefore, the last step of the protocol con- 
sists of applying the transformation [$fy ) —> ) 
described above in order to finalize the map \ty 



(k) 



l^^). It is important to remark that our protocol is fully 
deterministic. Despite the randomness of the photon- 
counter outcome, the determinism is recovered by choos- 



ing a different transformation \^f\, 



, |,T,( fc ) 



m/) for each 



I, such that the protocol always ends up in the final state 



I* 



(fc)\ 
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